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Abstract 

The Kravchuk and Meixner polynomials of discrete variable are introduced for the 
discrete models of the harmonic oscillator and hydrogen atom. Starting from Ro- 
f — , drigues formula we construct raising and lowering operators, commutation and an- 
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ticommutation relations. The physical properties of discrete models are figured out 
through the equivalence with the continuous models obtained by limit process. 
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^ ' 1 Introduction 

The method of finite difference is becoming more powerful in physics for differ- 
ent reasons: difference equations are more suitable to computational physics 
and numerical calculations; lattice gauge theories explore the physical prop- 
erties of the discrete models before the limit is taken; some modern theories 
have proposed physical models with discrete space and time [1]. 

In recent papers [2], [3] we have presented the mathematical properties of hy- 
pergeometric functions of continuous and discrete variable. We have worked 
out general formulas for the differential/difference equation, recurrence re- 
lations, raising and lowering operators, commutation and anticommutation 
relations. The starting point is the general properties of classical orthogonal 
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polynomials of continuous and discrete variable [4] of hypergeometric type, in 
particular the Rodrigues formula from which the raising and lowering opera- 
tors are derived. Similar results were obtained with more sofisticated method 
using the factorization of the hamiltonian [6] [7]. In those papers the term 
"oscillator" is used for the hamiltonian that has equally spaced eigenvalues. 

In this paper we study the physical properties of two simple examples, the har- 
monic oscillator and the hydrogen atom on the lattice. We make the Ansatz 
that the discrete model has the same elements (hamiltonian, eigenvalues, 
eigenvectors, expectation values, dispersion relations) as the continuous one, 
provided the difference cqiiation, the raising and lowering operator the com- 
mutation and anticommutation relations become in the limit the equivalent 
elements in the continuous case. To this scheme we can add the evolution 
of the fields with discrete time, using some difference equation that replace 
Heisenberg equation. [8] 

Some applications of discrete models have been presented elsewhere. Bijker 
et al. [9] have apply the SO (3) algebra of the Wigner functions to thcone- 
dimensional anharmonic (Morse) oscillator. Bank and Ismail have apply La- 
guerre functions to the attractiv Coulomb potential [10] 



2 The quantum harmonic oscillator of discrete variable 



We start from the orthogonal polynomials of a discrete variable, the Kravchuk 
polynomials K^^ (x) and the corresponding normalized Kravchuk functions [1] 

K(f\x)=d-'./^k^^\x), (1) 

where rf^ = n\{N-nyi P^)"' ^ normalization constant, p(x) = ^^lf^r-^(P9)" is 
the weight function, with p > 0, q > 0, p + q — 1, x — 0,1, . . . N + 1. 

The Kravchuk functions satisfy the orthonormality condition 

TV 

Y.K^\^)Ki?\^) = ^nn', (2) 
x=0 

and the following difference and recurrence equations [2]: 

^Jpq{N -x){x + l)K^\x + 1) 

+ ^pq{N -x + l)xK^^\x - 1) + [x{p -q)-Np^n\ K^^\x) = 0, (3) 



2 



^Jpq{N- n){x~^)K^l,{x) 

+ ^Jpq{N -n+l)nK^^l^{x) + [n{q -p) + Np-x\ K^J'\x) = 0, (4) 

From the properties of the Kravchuk polynomials we can construct raising and 
lowering operators for the Kravchuk functions [2] 

L+{x, n)K^^\x) = pq{x + n - N)K^\x) 

+ ^Jpq{N - X + l)xK^\x - 1) = ^Jpq{N - n)(n + l)K%{x), (5) 



L-(x, n)K^\x) = pq{x + n- N)K^p\x) 

+ ^pq{N-x){x + l)K^P\x + 1) = ^pq{N-n + l)nKi^'l,{x). (6) 

The raising operator satisfies 

= ^ ^'"^"i^^ "^' UL^i^^n- 1 - k)K^\x), 



where Kq (x) — y ^'\^l^_^y^ is the solution of the difference equation 

L-{x,0)Ki''\x)^0. 



It can be proved that the raising and lowering operators L'^{x, n) and L {x, n) 
are mutually adjoint with respect to the scalar product (2). 

If we define the difference equation (3) as the operator equation H {x,n)K^P\x) = 
, we can factorize this equation as follows 

L~^{x, n — l)L~{x, n) = pq{N — n + l)n + pq{x + n — 1 — N)H{x, n), (7) 
L~{x, n + l)L'^{x, n) = pq{N -n){n + 1) +pq{x + n + l- N)H{x, n). (8) 



Now we make connection between the Kravchuk function and the Wigner 
functions that appear in the generalized spherical functions [4] 

{-^y"-"'d^mmm = K^\x), (9) 

where j — N/2 , m — j — n, m! — j — x, p — sin^ (/9/2) , q — cos^ (/3/2) . 

Then formulas (3) to (8) can be written down in terms of the Wigner functions, 
namely, 
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- sin/3^(j + m'){j -m' + 



+ ^ sin (3^ {j - m') {j + m' + + (m - m' cos (3)d^^^^,{(3) = 0, 

(3a) 

^ sin + - m + 

+ ^ sin (j - m) (j + m + (/?) - {m! - mcos (5)d'^^^,{(5) = 0, 

(4a) 

L+(m',m) <^,(/3)=sin2 ^(m + m') <„,(/?) 
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+ - sin/3 J(j - m')(j + m' + 1) <„,+i(/3) 



\ sin /?V(i + ^)0'-^ + l) (/?) , 



(5a) 



L (m', m) (i^^^,(/3) =sin2 ^(^ + »^') 



+^ sin/?^(i + m')(i-m' + l) <^,_i(/?) 
= ^ sin/3^ (j - m)(i + m + 1) (i^^+i,^,(/?). 

(6a) 

Notice that (3a) and (4a) are equivalent if we interchange m <-> m' and take 
in account the general property of Wigner functions 

The same property is satisfied between (5a) and (6a). 

Exprcsions (5a) and (6a) can be obtained directly from the properties of 
Wigner functions. In fact, it is known [4, formula 5.1.19] that 

+ - V(j-m)(i + m + l)d;,^,,„,(/3), (11) 



12) 
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The last equation (11), after interchanging m m' and using (9), can be 
transformed into 



d ,i , n\ ~ cos /? „■ , ^, f—. — ; — , ^, 

^<m'(/?) sin/3 ^ ^.^'^^^ = V (.^ + )0 - ^' + 

(12a) 

Combining (11) and (12a) we obtain (6a) and by similar method we obtain 
(5a) In order to give a physical interpretation of the difference equation and 
raising and lowering operators for the Kravchuk functions we take the limit 
when N goes to infinity and the discrete variable x becomes continuous s . 

First of all, we take the hmit of Kravchuk functions. We write 



1 1 

' n\{N -n)\{NpqY\ [ NYq^-"" 1 ^2 

N\{pqy 2"(n!)2J \x\{N - x)\] 

1.. . , 1/. 



n 



Npq 



n 

n\k^\x) 



' I , ^ e-^'l EJs) = V'n(s), (13) 

where the last braket becomes the weight for the Hermite functions and the 
functions '0„(s) are the solution of the continuous harmonic oscillator (up to 

the constant (2A^pg) "^^4 . 

In order to get the continuous limit of (4) we multiply it by 2^^ ^J2Npq and 
substitute x = Np + \/2Npq s ; after simplification we get 



In the limit iV — > oo this equation becomes the familiar recurrence relation 
for the normalized Hermite functions 

^J2(^^^)^r^+l{s) + ^/Q^^r^-l{s)^2s^r^{s) (14) 

Before we take the limit of (5) we redifine the raising operator, substracting 
from it one half equation (3), namely. 
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L+{x,n)KjP\x) = ^!^[{x - Np) + n{p - q)]KlP\x) 

-^pq{N -x){x + l)KjP\x - 1) + ^Jpq{N -x + l)xK^P\x - 1) 
'pq{N - n){n + 1)kI^Ii{x) (15) 



We divide this expression by hy/2Npqy/Npq where hy/2Npq — 1. After sim- 
phfication we get 



-^L-^(x,n)K^\x)^^ 



s + 



n{p — q) 



1 

'2h 



A 



/ ^Nq'^Nqjy- 



I 2q 



s]Kn{x-l) 



n 



1--) in + l)K^Six) 



1 



71 ~ ^ ^'^'"^^ + /i) - A{s - /i)]| 



(16) 



Similarly from (6) we get 

^ L-(x,n) -L {, + I j ^„(,) = V5 (17) 

Therefore the raising and lowering operators for the Kravchuk functions be- 
come, in the limit, creation and annihilation operators for the normalized 
Hermite functions. 

We still have an other connection between the raising and lowering operators 
of Wigner functions with the generators of the SO (3) algebra. 

Prom (5a) and (6a) we define 
AUi^M ^ ^ LUl^,iP) = ^ U + m){j^-m + i ^ d?^_,,M. (18) 
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2j 



Multiplying both expressions by the spherical harmonics Yj^' , adding for m' 
and using the property of Wigner functions 



(20) 



we obtain 



A+v - U + m){j -m + 1) 



jm 



(21) 



A. 



(i - m){j + m + l] 



^ j,m+l 



1 



2j 



T V- 



where J+, J_ are the generators of SO (3) algebra. 
For the commutation relations of these operators we have 



n 



(22) 



(23) 



we substitute (20) in this expression and then take the limit: 



AA+ 



n 



„ , mm „ 



V'n(s) 



For the anticommutation relation we have from (7) and (8): 

[AA^ + A+A) Y^^ = ^ (j(j + 1) - m^) Y^^ =\{^- -^z) ^jrn. 
Again substituting (20) and taking the hmit 



(24) 



{AA+ + A+A)dl^,{P) 

= |(2n + 1) - y I di^,iP) — ^ (aa+ + a+a) M^) = (2n + 



This correspondence suggests that the operator algebra for the quantum har- 
monic oscillator on the lattice is expanded by the generators of the SO (3) 
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groups. The commutation relations [,/+, J_] = 2,]^ play the rol of the Heisen- 
berg algebra, and the anticommutation relation multiply by huj/2 play the rol 
of the Hamiltonian. In order to complete the picture we define the position 
and momemtum operators on the lattice as follows: 



X 



\] + m'){j - m' + I) j ij -m'){j + m' +1) 

2^ "m,m'-l \P) Y 2^ "m,m'+l \P ) 



with dispersion with respect to the state K^\x) 

(AXf^(x') ^^(aA+ + A+A) =_^f2n + l-^^ 
^ ^" \ /n 2Mcu \ /n 2Mco \ j\ 

from which the uncertainty relation follows: 



(AX)„(AP)„ = -(2n + l 



The eigenvalues of the Hamilton operator on the lattice are connected with 
the index m = j — n of the eigenvectors c^mm'(/^)- These eigenvalues are equally 
separated by hu but finite (m = — j, . . . + j). The eigenvalues of the position 
operator on the lattice are connected with the index m' = j — x of d-'^^,{l3). 

These eigenvalues are equally separated by \f^^ but finite {m! — —j, ...+]). 
Therefore the Planck constant h plays a role with respect to the discrete space 
coordinate similar to the discrete energy eigenvalues. 



3 Wave equation for the hidrogen atom with discrete variables 

Our model is based on the properties of generalized Laguerre polynomials as 
continuous limit of the Meixner polynomials of discrete variable. 

We start from the generalized Laguerre functions 

rn{s)^d-'4^)L-{s) (25) 
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with a > —1, (i^ ~ r{n + a + l)/n!, pi(s) = s"+-^e that satisfy the or- 
thonormahty condition 



(26) 



from the differencial equation and Rodrigues formula for the Laguerre poly- 
nomials [4] we deduce the following properties for the Laguerre functions (25) 



i) differential equation 

ii) Recurrence relations 



A _ 1 
s 4 



a 



i;:{s)^0, X^n + -{a + l). (27) 



-v/(n + a + l)(n + l)C+i(«) 

-^(n + a) nC-i(s) + (2n + a + 1 - s) ^/'^(s) = 0. (28) 

iii) Raising operator 

L+(s, n) = -^(2n + a + 1 - s) - C(«) 

= -^(^ + 1) + « + (29) 

iv) Lowering operator 



L-{s, n) C(^) = --(2n + « + 1 - 5) Cl^) + s- C(s) = 

= -Vri(ri + «)C-i(^) (30) 
from (29) and (30) we get the factorization of (27) 

L-{s,n + l)L+{s,n) = (n + l)(n + a + l) -s^ |^ + ^ - 1 - (31) 

L+(s, n - 1) L-(s, n) = n(n + a) - <^ — + - - (32) 

Notice that (27) corresponds to some self adjoint operator of Sturm-Lioville 
type. Also (29) and (30) are mutually adjoint operators with respect to the 
scalar product (26). 

Similarly, in the discrete case, we defined the normalized Meixner functions 

M^'^)(x) ^ <^v^ "^^'^H^) (33) 
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where m^''^^ {x) are the Meixner polynomials, 
2 n!r(n + 7) 



,2 ^ n!r(n + 7) . . ^ ^^r(x + 7 + 1) 

/."(i-/.)7r(7)' ^^^"^^ r(x + i)r(7) ^ 



and 7,/x are some constants 0</u<l, 7>0. The functions (33) satisfy 
the orthonormality condition 

00 -1 

X: M^^''^) (x) Mi7'^) (x) -^^^^ = (34) 

and the following properties: 
i) Difference equation 



/ Mx + 7)(x + l)(x + 7) 
X + 7 + 1 



+ -\//x(a: + 7)xM„(x - 1) - [^^{x + 7) + x - n(l - /x)] M„(x) = 0(35) 

ii) Recurrence relation 

- ^/x(n + 7)(n + l)M„+i(x) - ^//(n + 7 - l)nM„_i(x) 

+ (yux + //n + /U7 + n — x) Mn{x) — (36) 

iii) Raising operator 

L+(a;, n) M„(a;) = -//(x + 7 + n) M„(x) + ijl{x + ^)x M„(a; - 1) 

= ^/x(n + 7)(n + l) M„+i (x) (37) 

iii) Lowering operator 

L-{x, n) Mr,{x) = -^Ji{x + 7 + n) M„(x) 



M^ + 7)(x + l)(a; + 7) ^^^^ ^ 
V x + 7 + 1 



= -<^//(n + 7-l)nM„_i(a;) (38) 



A redefinition of (37) and (38) can be obtained substracting one half the 
difference equation (35) from both: 
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L~^{x, n)Mn{x) = — - (2n/u + n'j + n{l — n) + {fi — l)x) Mn{x) 



+i J ;M^M„(x - 1) - f'"^"j!,7;f ^"' m„(x + 1) } (37a) 

L-{x,n)Mn{x) = (2n// + /X7 + n(l - /x) + - M„(x) 



+5 ^M£±ll^:±iK^M„(, + 1) _ {^i^M^i,. - 1) \ (38a) 



It can be proved that the difference equation (35)corresponds to some self- 
adjoint operator of Sturm-Lioville type. Also the raising and lowering oper- 
ators (37) and (38) are mutually adjoint with respect to the scalar product 
(34). 

The anticommutation relations for the raising and lowering operators (37a) 
and (38a) are 



- |L"'"(x, n — 1)L (x, n) -|- L [x^n + \) {x^n)\ 



^(/x7+(// + l)n+(//-l)x)' 



A (x + 7)^(x + l)(a; + 7 + l)(x + 2) . 
[x + 7) (a; + 1) + y^(a; + 7)a;(a; + 7-l)(a;-l) (e^") 



yM^ + 7)xM„(x-l)i (39) 



The commutation relations read 

L+(x, n - 1) L"(a:, n) - L~{x, n^l) L+(a:, n) = -/x(2n 7) (40) 

In order to make connection between the Meixncr functions of discrete variable 
(33) and Laguerre functions of continuous variable (25) we substitute 7 — 
a-l-l, fj, — 1 — h, X — and then take the limit — > 0. 
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The explicite expression of (33) in terms oi a,h , and s reads 



^r(n + 7)^^' r(x + l) n\ 



r(a; + 7 + l)m(7'^)(x) 



1 - hY+^n\ 



^ r(n + a + 1) 
where we have used the asymptotic expansion 



n ^} 



r(x + 6) 



l + 0(- 



X ^ oo 



In the hmit we obtain (25), namely 



ni 



Also 



h-^o Vr(n + a + l) 

X—^00 ^ ' 

hx—*s 



(36) 
(37a) 
(38o) 

(39) 



(28) 
(29) 
(30) 



h^o 2 



o [(30) + (31)] 



In order to make application to the hydrogen atom we take the reduced radial 
equation [5] 



dp'' 



+ 



pi p^ 



u{p) = 



(41) 



where 



^J8M\E 



h V2|£;| 



If we compair (41) with (27) we must have p = s, + = (a^ — l)/4, hence, 



q; = 2/ + 1, and u — n + 



n + I + 1; since I is integer, a and v should 



be also integer. The energy eigenvalues are given by 

1 M{Ze'^f 1 



(42) 



For fixed u we still have degeneracy for / = 0, 1, ... — 1 
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The corresponding eignvcctors are given by (25), which after the substituion 
a — 21 + 1, and n — v — I — 1, become 



(^ + 0! 



) 



(43) 



Prom the connection between the Meixner and Laguerre functions given above, 
we can make the ansatz of a discrete model for the hydrogen atom where the 
reduced radial equation is substituted by the difference equation (35) with 
7 = q; + 1 = 2/ + 2, n = 1^-1-1, a; = 0, 1,... 

The raising and lowering opertors (29), (30), respect to the radial quantum 
number n , are substituted by the raising and lowering operators (37) and (38) 
with respect to n or = n + / + 1. 

The anticommutation relations (32) + (31), which is proportional to hamil- 
tonian of the hydrogen atom (for the radial part), is substituted by the an- 
ticommutation realtions (39). The commutation relations for the raising and 
lowering operators defining the Lie algebra of the SU(1,1) group, are substi- 
tuted by equation (40). 

The expectation value of the discrete variable x with respect to the Meixner 
functions Mj^'^^{x) is 



that can be calculated with the help of the recurrence relation (36). 

Finally the term /(/ + 1) in the Hamiltonian for the continuous and discrete 
case can be interpreted as the eigenvalues of the Casimir operator for the 



where the operators L+, L_ y L3 are acting on discrete space (see formulas 
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so (3) group. 




(5) (6)). 
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